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Detection  of  Targets  in 
Noise  and  Clutter 


1,  INTRODUCTION 

The  probability  density  function  (after  detection  and  video  integration)  for  a 
fluctuating  target  in  gaussian  noise  has  been  known  for  some  years  now.^  Similarly, 
results  have  also  been  obtained  for  the  case  of  a fluctuating  target  immersed  in  log- 

2 3 

normally  distributed  clutter.  ’ However,  it  does  not  appear  that  any  results  have 
been  obtained  for  the  case  of  a fluctuating  target  immersed  in  both  noise  and  clutter. 
Therefore,  in  this  paper  we  wilt  calculate  the  (post  square-law  detection)  probabil- 
ity density  function  for  the  sum  of  N pulses,  after  video  integration,  for  the  case  of 
a target  immersed  simultaneously  in  gaussian  noise  and  log-normally  distributed 
clutter.  We  shall  consider  principally  the  case  when  the  clutter  power  is  of  the 

same  order  of  magnitude  as  the  noise  power,  because  when  the  clutter  power  is 

2 

very  much  larger  than  the  noise  power  the  results  in  are  appropriate,  whereas 
when  the  noise  power  greatly  exceeds  the  clutter  power  the  results  in^  are  adequate. 


(Received  for  publication  21  July  1977) 

1.  Marcum,  J. , and  Swerling,  P.  (1960)  Studies  of  target  detection  by  pulsed 

radar,  IRE  Trans,  on  Inform  Theory  IT -6:2. 

2.  Trunk,  G. , and  George,  ,S.  (1970)  Detection  of  targets  in  non-gaussian  sea 

clutter,  IEEE  Trans,  on  Aerospace  and  Electronic  Systems  AES-G:620. 

3.  Schleher,  D.  (1975)  Radar  detection  in  log-normal  clutter,  IEEE  International 

Radar  Conference.  IEEE  Press  (Publication  No.  75  CH0938-1  AES) 

New  York. 


Wo  shall  present  detailed  numerical  results  for  the  probability  of  a false 
alarm  and  the  probability  of  detection  for  both  Sworling-l  and  Swerling-J  ti;rgets. 
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2.  TIIKORKTK'.M,  ANALYSIS 


2.1  Tarjict  and  tilutler  tion.siant  1 Dim  I’uIsf-lo  PuIsc  but  lurying 
From  Scun-lo-Sciin  (Swcrliiig  1) 

Let  U.S  suppose  we  have  a target  immersed  in  gaussian  noise  and  log-nornially 
distributed  clutter,  and  this  target  is  to  bo  detected  using  an  N-pulse  radar  burst. 
If  we  as.sume  thac  both  the  target  and  the  clutter  are  constant  from  pulse-to-pulse 
but  vai-y  from  scan-to-scan,  then  the  single  pulse  return  Re'^  can  be  written  as 


Re 


ip 


p + re 


(1> 


£ I 
» 
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where  p represents  the  target -plu.s-clutter  and  r expliili)  is  a complex  phasor  rep- 
resenting the  noise,  which  is  as.sumed  to  vary  from  pulse-to-pulse,  and  have  the 
probability  density  (for  r 5 0) 


P_(t',  <t>)  P exp 


1;^) 


1.2) 


for  mathematical  convenience,  we  assu-.ne  that  the  signal  is  to  be  detected 

using  a square-law  device.  We  are  therefore  interested  in  the  probability  dislribu- 
2 

tion  of  the  square,  R , of  the  envelope  of  the  received  signal.  If  we  define 
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p(ylx)  = exp(-x  ••  y)  Iq  (2-J^  for  y > 0 
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= 0 for  y < 0 , (3) 

where  ) is  the  mudiried  Bessel  I'unctLon.  The  characteristic  function  Cjisjx) 

associated  with  p(ylx)  is 

oc 

r 

C^islxlsJ  dye'^^p(yix) 
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00 

- e*’' Jdy  (4) 

0 

= <s  -t  t)"^  exp  ^-x  + 

For  N-pulse  detection,  we  detect  the  sum  of  N different  values  y^  of  the  signal 
plus  noise  vttriahlos.  That  is, 


V ^ yj  ^ ^2  ^ y.i  • ■ • " y.N 


(5) 


If  the  noise  samples  are  statistically  independent,  then  the  characteristic  function 
of  the  sum  of  the  N samples  is  the  product  of  iheir  individual  characteristic  func- 
tions, Ttierefore,  tlie  characteristic  function  Cj^(s',  x)  associated  with  T is 


Cj^^islx)  = (s  + exp 


^-Nx  + 


Nx 

s 1 


(G) 


Now  because  the  target  and  clutter  fluctuate  from  scan-to-scan,  we  must 

*1  2 

therefore  average  (fi)  over  the  variable  x - P*  • That  is 


00 

Cj^(s)  - j dx  p(x)  Cj^(s|x) 

-•X) 


* dxp(x)  exp(-Ny.  + 

L 


(7) 


where  p(x)  is  the  probability  density  of  the  random  variable  x.  We  shall  calculate 
this  later.  Once  C^(c)  is  known,  the  probability  density  of  Y = y s y2  • • • is 


Pn<Y) 


/ 


ds  Cjg(s)  e' 


= J dx  p(x)  e Gj^,(x,  Y)  , 


Gn(x.y)=^.-j  'i^(2.my)-.  (^1 


y'(N-1)/2 


(2yNxY) 


The  result  in  Eq.  (9)  follows  from  the  (N  - 1)  fold  differentiation  of  both  sides  of 
the  known  result 


"Ai  / 


ds  exp  (sY  4 ^ 


In  order  to  evaluate  Eq.  (8)  for  p(y),  we  must  first  obtain  p(x).  We  recall 
that  p is  the  sum  of  the  target  and  clutter  and  can  be  written  as 

p = p exp  (ix)  + n exp  (itf/)  (11 

where  p exp  (iX)  is  a phasor  representing  the  target  return  and  n exp  (ii/)  is  a 
phasor  representing  the  clutter  signal.  For  a Rayleigh  target,  we  know  that 


/ 2 \ 

Pj<w,  X)  = exp  j - j for  p i 0 

2^0t  2aj  y 


= 0 for  p < 0 


and  the  log-normal  clutter  is  governed  by  the  probability  density 


= (2ir)'^^^  (on)*^  exp 
= 0 for  ri  < 0 

The  quantity  r)  In  (13)  is  the  median  value  of  t|.  Also,  the  standard  deviation  a 
° 2 3 

has  been  estimated  ' to  be  0.7  for  low-se^-state  clutter,  1.  1 for  ground  clutter, 
and  as  high  as  1. 6 for  hi(^h-sea-state  clutter. 

The  probability  density  for  the  envelope  of  a Rayleigh  plus  log-normal  phasor 

4 

is  given  in  Chapter  4 of  Beckmann  and  can  be  written  as 

dr  F(T)  Ij,  (S  for  p 2 0 (14; 

= 0 for  p < 0 , 

whe  re 


F(t)  = (2ir)  (oo^t)  ^ exp 


'"It) 


2a‘  2a 


2 2 

If  we  recall  that  x = p /^°n’  easy  to  show  that 


Ptc<*>  = «n  e‘‘”‘|dTF(T) 


for  X 2 0 


where 


b ■ 


(15) 


We  now  use  (15)  in  (8)  and  set  n " and  After  some  manipulation, 

we  obtain 


p (y)  = 7e-Y(Yf  J ° 


= 0 for  Y < 0 


■i^(aBr)h 

8a^ 


= noise  to  clutter  ratio 


signal  to  noise  ratio 


Equation  (l.fi)  is  the  probability  density  for  the  sum  of  N pulses  (following  square- 
law  detection)  received  from  a slowly  varying  target  immersed  in  slowly  varying 
clutter  but  rapidly  varying  noise.  This  can  be  thought  of  as  a Swerling-1  target 
immersed  in  Rayleigh  noise  and  Swerling-l  clutter. 

As  a check  on  the  validity  of  (16),  we  can  study  the  limit  as  the  clutter  power 
vanishes  (jl -*  «).  In  this  case,  (17)  should  reduce  to  the  usual  Swerling-l  result. 
When  a it  is  clear  that  the  only  contribution  to  the  integration  over  n in  (17) 
comes  from  n very  near  iero.  Therefore,  we  can  expand  the  integrand  in  a 
Taylor  scries  about  rj  = 0.  This  gives 


^ f dn  [ I n^  (o/3r7^)1 

/3  -»  00  ^ L 8o^  J 


If  we  set  M(5)  = 1 in  (16)  we  readily  find  using  p. 


187  of  Bateman 


5 


that  for  Y £ 0 


p (Y) 

N li  ^co 


(N  - 2): 


(21) 


5 

where  7(N,  x)  is  the  incomplete  gamma  function  defined  on  p.  387  of  Bateman. 

Equation  (21)  is  precisely  the  well-known  result  for  a Swerling-1  target  in  Rayleigh 

0 

noise,  as  is  seen  by  comparison  with  the  results  in  Section  3,  5 of  Berkowitz. 

An  additional  check  is  obtained  by  studying  the  result  for  N = 1.  For  y £ 0, 

(16)  becomes  in  this  case 


Pl(Y)  = 


7 Y 

7 exp 

00  1 
f ±1  T 

n(27Y)^^^ 

1 

exp  ' 
1 

1 ^ 

_21 

. 2,  „ 2. 

J n (cfSn  ) 

(2jr)^/^  0(1  + y) 

J ^ 1 

Q 

1+7 

1 2(1  + 7)  ■ 

1 

8o^ 

(22) 


which  is  the  well-known  result  for  the  sum  of  two  Rayleigh  and  one  log-normal 

• random  variable.s.  Equation  (22)  is  also  appropriate  for  the  N-pulse  case  in  the 
limit  when  the  noise  samples  are  completely  correlated  from  pulse-to-pulse,  as 
opposed  to  the  case  we  have  considered  here  with  the  noise  samples  statistically 

• independent  from  pulse-to-pulse. 

1 

( 

! 2.2  Clutter  Varying  1 roin  Scan-to-Scan  but  Target  V arying  From 

i Pulse-to-Pulse  (Swerling-2  Target) 

i When  the  target  varies  so  rapidly  that  its  cross  section  changes  from  pulse-to- 

pulse  (Swerling-2  Target)  but  the  clutter  is  constant  from  pulse-to-pulse  (where 
i the  clutter  is  assumed  to  vary  from  scan-to-scan)  we  can  replace  Eq.  (1)  by 

V 

j R e*'’  - u + V e*®  (23) 

I where  u represents  the  clutter  and  v exp  (i6)  is  a phasor  representing  the  target 

j plus  noise.  For  a Rayleigh  target  immersed  in  Rayleigh  noise  (that  is,  the  noise 

I envelope  is  Rayleigh,  although  the  quadrature  components  are  gauss  ian),  the  joint 

I probability  density  of  v and  6 is  given  by  (for  v £ 0) 

f 


5.  Erdelyi,  A.,  et  al  (19.64)  Tables  of  Integral  Transform  (Bateman  Manuscript 

Project),  McGraw-Hill,  New  York. 

6.  Berkowitz,  H.  (1965)  Modern  Radar.  Wiley,  New  York. 
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Pjn(v,6)  = 


0/2  2v 

27r(On  ' 


exp 


(24) 


2 2 

where  is  the  noise  variance  and  is  the  target  variance.  VVe  are  again 
interested  in  obtaining  the  probability  density  of  the  square  of  the  envelope  of  the 
total  received  signal.  Let  us  define 


2 2 


y 


R 


2(o2  . a2) 


y 

1 + a 


where  y is  the  same  variable  defitied  following  Eq.  (2)  and  u is  defined  in  (20).  If 
we  temporarily  treat  u as  a constant  because  it  does  not  vary  from  pulse-to-pulse 

4 

(but  only  from  scan-to-scan),  it  is  straightforward  to  show  that 


p(y' 


I z)  = exp  (-Z  - y')  Iq  |^2(2y')*^^j  for  y’  s 0 


= 0 for  y'  < 0 


(25) 


which  gives  a single-pulse  characteristic  hmction 

Cj^(sl  z)  = (s  + 1)'^  exp  ^-z  + g- ^ ^ ^ . (26) 

For  N-pulse  detection,  we  detect  the  sum  of  N different  values  yj^.  In  partic- 
ular 


Y'  ^ y'^+y^  + ^yj^,  . (27) 

If  we  assume  that  the  target  plus  noise  variables  are  statistically  independent  from 
pulse-to-pulse,  the  N-pulse  characteristic  function  is  the  product  of  the  single 
pulse  ones  and  we  get 


I 


i 


16 


(28) 


Cj^j(s|  z)  = (s  + 1)  ^ exp  1^"  Nz  + 

Because  the  clutter  varies  from  scan-to-scan,  we  must  next  average  (28)  over 
the  clutter  probability  density.  This  gives 


Cj^(s) 


W 

I 


dz  p^(z)  exp  (-Nz 


(s  + 1) 


N 


(29) 


where  p (z)  is  the  probability  density  of  the  square  of  the  clutter  envelope.  If  we 

write  the  clutter  u as  the  random  phasor  u = n exp  (h^;),  it  is  evident  that  p in,  il)  is 

2 2 2 

given  by  (13).  In  terms  of  the  variable  z = u /2(o^  + e^),  this  becomes 

p^(z)  = ^2(2)r)^^^  czj’^  exp 

- 0 for  z < 0 
2 2 2 

where  h = (o^  + If  we  use  (30)  in  (29),  we  get 

Nz  + -ii£_  -ij£i2h.Sl 


C^.,(s)  = 


2(27r)‘^^  ct(s 


/dz 


exp 


tn^(2hz) 


6a 


for  z a 0 


(30) 


so  that 


ii+ioo 


PN^^^"27ri  j ® Cj^(s) 


A-loo 


is  given  by 


•Nz  - 


fn  (2h  z) 
8a^ 


Gj^(z,  Y')  for  Y'aO 


(32) 


= 0 for  Y'  < 0 


f! 


where  Gj^(z,  Y')  is  given  by  (9).  Finally,  if  we  let  ? = Nz  we  can  rewrite  (32)  as 


« 2(2,)‘'2,  J ,5fN*iy2  ‘N-l[2"n  J 

o 


= 0 for  Y'  < G 


exp 


- € - 


80“ 


for  Y'  > 0 


(33) 


As  a check  on  the  validity  of  (33),  we  can  consider  the  limiting  case  when  the 
clutter  is  absent  (h  ->■  oe  ).  In  this  case,  it  is  evident  that  the  dominant  contribution 
to  the  5 integral  comes  from  5 near  zero.  Therefore,  we  may  expand  the  integrand 
in  a Taylor  series  about  C = 0-  If  this  is  done,  the  5 integral  may  be  performed 
and  we  find  (for  Y'  2:  0) 


Pm(Y')  = 


(Y')N-I  e-^' 

(N  - D! 


which  is  identical  with  the  well-known'’  result  for  the  N-pulse  sum  for  a Swerling-2 
target  in  Rayleigh  noise. 


3.  PROBABILITY  OK  FALSE  ALARM  ANO  OETECTION 

The  probability  of  a false  alarm  is  the  probability  that  the  received  signal 
exceeds  some  threshold  level  Y when  there  is  no  target  present.  In  the  limit 
when  the  target  is  absent,  it  is  readily  demonstrated  that  for  Y 2:  0,  (16)  becomes 


^N-1  exp 


-? 


‘AW) 


8a‘^ 


(3  4) 
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Therefore,  the  probability  of  a false  alarm  for  an  N -pulse  detection  system  is 


'■i-- 


2(27r)‘-‘"a 


d|  If^.i  [2(Y?)^/^J  exp 


- P . 


(?4i) 


(?)' 


>(N+l)/2 


8a 


(35) 


The  probability  of  detecting  s Swerling-l  Target  in  noise  and  slowly  varying 
ter  can  be  obtain 
is  present.  This  is 


clutter  can  be  obtained  from  (16)  and  is  the  probability  that  Y — when  a target 


f -Y  /N-l)/2  f d?  I [2(yc)^^^]  M(C)e'‘^+^>^ 

Pd,  = y / dY  e ^ (Yf^  J (36) 

Y o 

o 


Also,  the  probability  of  detecting  a SwerUng-2  target  in  noise  plus  slowly  vary 
:lutter  is  the  probabi 
fore,  from  (33)  we  have 


ing  clutter  is  the  probability  that  Y'  & Y^/(l  + u)  because  Y'  = Y/d  + a).  There- 


d2 


-W-  f dYe-Y(Yi^-l)/2  f 

2(27r)  a J 


Is^ 

I'to 


(37) 


6t  1. 


It  is  easy  to  demonstrate,  that  P..  = P for  N = 1. 

d 1 da 


In  the  next  section,  we  will  present  numerk  results  for  P,,  P^,  and  P^2' 
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1.  MMKRICAL  RKSliLTS 


Before  presenting  our  results  it  is  useful  to  reemphasize  once  again  our 
definitions  for  signal-to-noise  and  noise-to-clutter  ratios,  since  they  may  differ 
somewhat  from  those  used  hy  other  authors.  They  are 
^2 

a = — signal-to-noise  ratio, 

2 

a 

,3  = — ^ = noise  to  median-clutter  ratio, 

N = number  of  independent  pulses  integrated, 

Pj.  - probability  of  a false  alarm, 

P = probability  of  detecting  a Swerling-l  target  in  noise  varying  from 
pulse-to-pulse  and  clutter  varying  from  scan-to-scan, 

P = probability  of  detecting  a Swerling-2  target  in  noise  varying  from 
pulsc-to-pulse  and  clutter  varying  from  sc.in-to-scan. 


4.1  I'ulse-  VlHrm  Probubilit) 

By  using  (3  5),  we  have  calculated  the  probability  P.  of  a false  alarm  as  a 

2 2 ^ 

function  of  the  threshold  Y = R /2o  . for  a number  of  different  values  of  the 

o o'  n 

noise-to-clutter  ratio  3,  These  results,  with  the  standard  deviation  a of  the  log 
normal  distribution  set  equal  to  0.7,  are  presented  in  Figures  1-4.  We  observe 
that  the  probability  of  a false  alarm  is  significantly  affected  by  the  value  of  3.  It 
is  even  mare  strongly  affected  by  the  value  of  a as  can  be  seen  from  the  resu  » in 
Figures  5 and  0. 

From  the  nature  of  the  curves  in  Figures  1-6  it  is  clear  that  for  large  values 
of  (that  is,  small  false  alarm  probabilities),  P^  is  a function  of  (lY^/N.  It  is 
demonstrated  in  Appendix  B that  for  Y^  large,  we  may  approximate  (35)  by 


— erfc 
2 


,,3/2 
2 ' c 


I ^ a 


(38) 


where  erfc  is  the  complementary  error  function  defined  by 


erfciz) 


2 


(.)^/2 


/ 2. 
exp  (-Z  ) 


(39) 
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Pf 

(exact) 

(calculated  from 
Equation  (38)) 

2.  5a  X lO'^’ 

2.  57 'X  10’'' 

4.  20  X 10"^ 

4.  20  X 10‘^ 

1. 75  X 10'^’ 

1.  79  X lO”'’ 

4.  10  X 10'^ 

4.  20  X lO*”^ 

2.  0 X 10'^ 

1 

2.  12  X lO'^ 

t.2  l*riiliabilil\  olDctfrlion  for  S»erlin(;-2  larpol.s 

We  now  present  some  numerical  ■'esuUs  for  the  probability  of  detecting  a 
Swcrltng-2  target  immersed  in  Rayleigh  noise  which  varies  from  pulse-to-pulse 
and  log-normally  distributed  clutter  which  fluctuates  from  scan-to-scan.  By  using 

Eq.  (37),  we  have  calculated  the  pi’obability  of  detection  P for  false-alarm  prob- 
-4  -6 

abilities  of  10  and  10  . These  calculations  have  been  performed,  assuming 

0-0,7,  for  a number  of  values  of  ;3  and  N,  and  are  presented  in  Figures  7-lG. 

For  convenience  we  have  also  presented  in  Table  2 the  values  of  the  signal -to-noise 
ratio  <>  required  for  ^ ® number  of  different  values 

of  noise-to-clutter  ratio  and  pulses  integrated. 

The  effect  on  P^2  changing  the  standard  deviation  o of  the  log-normal  dis- 
tribution IS  illustrated  in  Figures  17  and  18.  We  note  that  a change  in  a from  0.  7 
to  1.  1.  results  in  more  than  one  order  of  magnitude  change  in  the  value  of  signal 
to  noise  ratio  a required  for  a given  probability  of  detection  P^2'  case  of 

vanishing  clutter  (3  -»  «>,  our  numerical  data  compares  well  with  the  data  in  refer- 
ence 1 for  a Swerling-2  target  immersed  in  Rayleigh  noise.  In  the  opposite  limit 

7 

of  vanishing  noise  |3  .*  0 we  can  compare  our  results  with  those  of  Trunk  only  for 
N = 1,  because  Trunk  has  considered  tlie  case  v/hen  the  clutter  varies  from  pulse- 
to-pulse  and  we  have  considered  the  limit  when  the  clutter  varies  only  from  scan- 
to-scan. 

From  Figures  7-lB,  we  note  that  large  signal-to-noise  ratios  are  usually 
required  for  high  detection  probabilities.  This  suggests  that  it  might  be  useful  to 
have  a large  n approximation  to  (3  7).  For  large  a,  it  is  shown  in  Appendix  C that 
(37)  can  be  approximated  as 


m = 0 


A comparison  in  Table  3 of  numerical  results  calculated  from  (40)  with  the  com- 
puted data  in  Figures  7-lH  indicates  that  (40  is  quite  an  excellent  approximation 
for  P^2  provided  oil/N  » 1.  This  encompasses  nearly  the  entire  range  of  practical 
interest  (that  is.  P^g  > 0.  .t).  Therefore,  in  quite  a number  of  practical  situations 
it  is  not  neces.sary  to  numerically  evaluate  the  integrals  in  (37),  and  Eq.  (40)  can 
be  used. 


7. 


T 


runk,  G.  (1971)  Further'  results  on  the  detection  of  targets  in  non-gaussian 
sea  clutter.  IEEE  Trans,  on  Aerospace  and  ITectronic  .Systems  AES-7:553. 
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Figure  11.  Probability  of  Detection 
(Sterling  2)  for  P,  = lO'"^.  p = 100, 
a --  0.7 


Figure  13.  Probability  of  Detection 
(Swerling  2)  for  P.  “ lO**’,  ;(  - 1, 

CT  = 0.  7 ' 


/3  ■ 0 316 


N • 3 


Figure  12.  Probability  of  Detection 
(.Sweri.ng  2)  for  Pf  = lO"**,  ;i  = 0,310, 
a = 0.  7 


Figure  14.  Prcibability  of  Detection 
(Swerling  2)  for  P,  - 10“^’,  ,3  = 3.10, 
a - 0.  7 ' 
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Table  3.  Comparison  of  Pd2  Calculated  From  Equation  (40)  With  the  Exact 
Result  for  a = 0.  7 and  Pj.  = 10“*> 


a 

N 

^o 

’^d2 

(exact) 

*^d2 

(calculated  from 

Equation  (40) 

1 

1000 

1 

390 

0.  67 

0.  68 

1 

1000 

3 

1170 

0.  89 

0.  89 

1 

1000 

r. 

23  50 

0.  968 

0.967 

1 

1000 

10 

3900 

0.  994 

0.  9932 

10 

100 

1 

44 

0.  65 

0.  65 

10 

100 

3 

123 

0.  875 

0.  876 

10 

100 

240 

0.  965 

0.  966 

10 

100 

10 

405 

0. 9922 

0. 9918 

1.3  Prolmbili(\  of  llolerlion  for  Swt*rltnp-I  Targols 

We  now  calculate  the  probability  of  detecting  a S'werling-l  target  immersed  in 
log-normal  clutter  which  varies  from  scan-to-scan  and  Rayleigh  noise  which  varies 

from  pulse-to-pulse.  By  using  (3f>),  we  have  computed  the  probability  of  detection 

-4  -0 

for  false-alarm  probabilities  of  10  and  10  . These  calculations  for  o - 0,  7 

and  a number  of  different  values  of  /3  and  N are  presented  in  Figures  19-28.  From 
Figures  19-28,  we  observe  that  large  signal -to-noise  ratios  o are  usually  reqxiired 
for  high  probabilities  of  detection.  For  large  o,  it  is  readily  demonstrated  that 


(41) 


A comparison  of  numerical  results  calculated  from  (41)  with  the  computed  data  in 
Figures  19-28  indicates  that  (41)  is  an  excellent  approximation  for  whenever 


Tlie  proof  parallel.s  that  in  Appendix  C.  In  particular,  observe  that  for  large  o 
the  integrand  in  (17)  can  be  expanded  in  a Taylor  series  about  rj  = 0.  This  gives 
M(f)  ■ 1.  Upon  setting  M(C)  = 1 in  (3G),  we  can  readily  do  the  f.  and  y 
integrations. 
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Figure  27.  Probability  of  Detection 
(Swerling  1)  for  P,  - 10'^,  (3  = 10, 
a --  0.  7 


Figure  28.  Probability  of  Detection 
ISwerling  1)  for  P,  = 10**^,  j3  - 100, 
o = 0,  7 
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It  is  interesting  to  note  from  Figures  19-28  that  unless  ;'i  > 10,  P^^  is  essen- 
tially independent  of  N.  This  is  quite  evident  from  physical  reasoning  because  for 
small  noise-to-clutter  ratios  neither  the  target  nor  the  clutter-plus-noise  fluctuate 
from  pulse-to-pulse.  Therefore,  summing  N-pulses  cannot  offer  any  advantage. 
This  is  also  evident  from  (41)  and  (88).  P rom  (38)  we  have  t.iat  for  large 


y 

o 


H(Pj.)N 


where  If  is  a function  of  Pp  If  we  use  this  result  in  (41)  we  obtain 


dl 


H(P^) 

Ojl 


(42) 


(48) 


which  is  clearly  independent  of  .N'.  For  large  values  of  the  results  in  (42)  and 
(43)  are  no  longer  appropriate  because  then  is  no  longer  large,  as  is  required 
for  (42)  to  be  valid. 


5.  DISCI  SSION 


The  results  we  have  presented  in  Sections  2-4  are  strictly  valid  for  a receiver 
which  consists  of  a square-law  detector  and  integrator.  However,  most  systems 
operating  in  an  environment  where  clutter  is  a consideration  are  likely  to  employ 
an  MTI  (moving  target  indicator)  to  reduce  clutter.  V/hen  an  MTI  is  used  our 
results  are  only  approximately  valid,  because  the  noise  samples  are  then  not 
statistically  independent  from  pulse-to-pulse,  as  we  have  assumed,  (This  state- 
ment applies  to  the  target -plus -noise  for  a Swerling-2  target.)  In  order  to  see 
this,  we  consider  a single  stage  MTI.  Let  n._j  be  the  noise  present  at  the  input  to 
the  MTI  at  time  (i  - 1)t,  n.  be  the  noise  at  time  ir,  and  n.^j  be  the  input  noise  at 
time  (i  + 1)t.  Then  the  MTI  noise  output  n^  at  time  iT  is  n.  = n._j  -n^  and  the  output 
at  time  (i  + 1)t  is  n.^j  = n.  - Therefore,  even  though  n.  _j,  n.  and  n.^j 

are  statistically  independent,  and  n.^^  are  not.  because  (n.  " "j) 

(n^  - n-^j))  = - (n^)  * 0.  Consequently,  the  noise  samples  at  the  output  of  the  MTI 
are  partially  correlated  from  pulse-to-pulse.  Thus  if  N of  these  pulses  are 
integrated,  the  correct  probability  of  detection  will  be  somewhere  between  our 
answers**  in  Figures  7-28  for  N = I (complete  correlation  of  all  pulses)  and 
N = (complete  decorrelation).  For  a single  stage  MTI,  the  pulse-to-pulse 
correlation  coefficient  is  1/2,  for  a two  stage  MTI  the  correlation  coefficient  is 
2/3,  etc.  This  suggests  that  the  N = I results  in  Figures  7-28  will  be  a reasonable 
approximation  when  a multistage  MTI  is  used,  even  if  N -pulses  are  integrated 
following  the  MTI. 


*Note  that  this  will  not  be  a problem  with  the  clutter  samples  because  the  clutter 
samples  have  been  assumed  to  be  completely  correlated  from  pulse-to-pulse, 
and  the  MTI  does  not  change  this.  The  clutter  stati.stics  from  scan-to-scan  are 
of  course  still  log-normal. 

**The  value  of  the  clutter  level  which  would  be  used  in  calculating  0 is  the  residual 
clutter  level  present  at  the  MTI  output  terminals. 
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Appendix  A 

Probability  Dansity  for  a Constant  Target 


Here  we  give  the  probability  denasfy  for  the  envelope  squared  of  the  sum  of 
N-independcnt  .samples  from  a constant  target  imniersed  in  Rayleigh  noise  which 
varies  from  pulse-to-pulse  and  log-normal  clutter  which  varies  fi-om  scan-to-scan. 
For  Y a 0,  the  result  is 


P(y)  = ■—  ,,  /g-  ■ r d?  e"'  G-.(E,  Y)  f 
{2iry'^o  J ^ J 


dii  exp 


K 


•<  /?  1 'V  2 -1 

jj(2£  - 2'^'  ^ a'^)N  ^ 

8a2 


2|  - a?l'2  + a2 


(.Al) 


1 

where  a - N “C/a  , C is  the  (constant)  target  voltage  and  (3,  Y,  Gj^  etc.  are  as 
defined  previously. 


Appendix  B 

Dtrivotion  o(  Equation  (38) 


If  we  let  n = Y/V  and  f = Y t,  we  can  rewrite  (35)  as 
' o o 


■r  Y^/  <"-/ 


exp 


For  Y » 1,  we  can  use  the  asymptotic  expansion  for  L ,(Z),  If  this  is  done  and 
o 2 2 ™ ^ 

we  also  make  the  substitution  t = and  n = v , we  can  rewrite  (Bl)  as 


Pr  = 


{2Y  ) 
a_ 


1/2 


f 2 jra 


I dv  J ^ exp  -Y^(x  - v)2  - ^ f n 

1 o L 


For  Y , 00  it  is  clear  that  the  integrand  of  the  x integral  is  sharply  peaked  about 

o 2 

X = V.  because  of  the  term  exp  (-Y„<x  - v)  ) . Therefore,  we  can  set  x = v in  the 

® 2 

entire  x integrand,  except  for  exp  [-Y  (x  - v)  ] to  get 
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(2Y  )^/2  op 


Pf- 


27ro 


J ^ 


exp 


In 


2o‘ 


« • £ 
-Y  z 


dz  e 


(B3) 


where  in  the  last  integral  in  (B3)  we  have  replaced  x - v by  z and  extended  the 
range  of  integration  on  z to  the  entire  real  line.  If  we  perform  the  z integral  and 
then  let  u = In  [v(2(lY  we  get 


\ f 


du 


(B4) 


where  p = (tla^)  In  (2;iY^/N).  Equation  (38)  follows  immediately  from  (K4), 
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Appendix  C 

Derivation  of  Equation  (40) 


For  lar-ge  a,  it  i.-s  convenient  to  rewrite  (37)  art 


Pd2  ^ - 


^i:^  / 


o o 


As  i»  - <jc , it  is  clear  that  most  of  the  contribution  to  the  integral  comes  from  t 
near  zero.  Therefore,  we  may  expand  in  a Taylor  series  about  f - 0 to 


obtain 


P,2  - 1 - 


1 r * dY  e'^  f r 1 2/20 3 ^'\1 

^ J ^N^rnr—  J T -;;2  • 

r\  n •— 


If  we  recall  that 


‘ 


we  see  that  (C2)  becomes 


l+< 

-/ 


dY  e'"^ 

(N  - 1)! 


Finally,  if  we  use  form  5 on  p.  134  of  Hef.  I,  we  can  evaluate  the  integral  in  (C3) 
to  yield  (40). 


Erdelyi,  A. 
Project), 


, et  al  (1954)  Tables  of  Integral  Transform  (Bateman  Manuscript 
McGraw-Hill,  New  York. 
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SI  Symbol 

Fannula 

length 

metre 

m 

meH 

kilogram 

kg 

timt 

eecond 

a 

tltctric  cufTvnt 

ampere 

A 

theimodynamlc  Itmptnture 

kelvin 

K 

•nouQi  of  tubitanco 

mole 

mol 

lumiaooi  ItttoMity 

cendele 

SUmBMEhTTAKV  UNITS: 

plene  ooglo 

redlan 

red 

•olid  uiilo 

elemdlan 

er 

DERIVED  UNITS; 

Aocrtontioa 

metre  per  eecond  equered 

... 

aWa 

•cdvlty  (of  • rodlooctlvo  oourc*) 

dUinttgrtUon  per  eecond 

... 

(dialategNtlonVa 

tnguUr  •ootUrtlion 

radian  per  eeoond  equarad 

... 

tad/a 

•ngulir  velocity 

radian  per  eecond 

... 

tad/a 

Ofee 

equere  m^ie 

... 

m 

denelly 

kllognm  per  cubic  metre 
fared 

... 

kgftn 

electric  cepecltence 

F 

A-aA^ 

electrical  conducting 

eiemene 

S 

AAf 

electric  field  rirength 

volt  per  metre 

V/m 

electric  Inductance 

betiiy 

H 

V‘efA 

electric  potential  difference 

volt 

V 

W/A 

electric  lealetance 

ohm 

V/A 

electromotive  force 

volt 

V 

W/A 

•neffy 

ioule 

1 

N*m 

entropy 

loule  per  kelvin 

force 

newton 

N 

kgaUB 

frequency 

heiis 

Ht 

(cydaVi 

illumlnenoe 

lux 

Ix 

Im/ro 

lumlnenco 

candela  per  aquara  matia 

odftn 

lumlnou*  flux 

lumen 

im 

cd^ 

magnetic  field  itrength 

ampere  per  metre 

... 

A/m 

magnetic  flux 

vreber 

Wb 

V>a 

magnetic  flux  denelty 

teale 

T 

WVm 

magnetomotive  force 

ampere 

A 

power 

Witt 

W 

preeaure 

paecal 

Pe 

N/m 

quantity  of  electricity 

coulomb 

C 

A*a 

quantity  of  heet 

loule 

1 

N>m 

latent  Intenaity 

watt  par  eteredien 

... 

W/ar 

epecinc  beet 

loule  per  kilogrem-kelvfn 

Pkg-K 

etreae 

paecal 

Pe 

Nhn  ' 

thermal  conducUvity 

watt  per  metre^k^vln 

... 

W/m-K 

velocity 

metre  per  lecood 

m/a 

viiGOilty.  dynamic 

paecal'eecond 

... 

Pa-a 

vlicMlty,  kinematic 

aquart  metre  pw  eecond 

... 

m/a 

voltage 

volt 

V 

W/A 

volume 

cubic  metre 

m 

wavenun^MT 

work 

reciprocal  metre 
loule 

i 

(weveym 

N>m 

SIPUnXES: 


Multiplication  Fedora 

Ptefix 

SI  Syq 

1000  000  000000-  XO** 

ten 

T 

1000  000  000-  10* 

filf* 

C 

1 000  000-  10* 

mage 

M 

1 000-  10> 

kilo 

k 

100  - 10> 

becto* 

b 

10  - 10* 

deka* 

da 

0.1  - 10-* 

decl* 

d 

0.01  - to-> 

Dentr 

c 

0.001  - t0*“> 

mllli 

m 

0.000  001  - 10** 

micro 

M 

0.000  000  001  « 10** 

nann 

n 

0.000000  000001  - 10*** 

ptoe 

iamlo 

P 

0.000  000  000  000  001  - 10*** 

f 

0.000000  000000  000  001  « 10*** 

atto 

a 

* To  be  ivolded  when  poMible. 


